🎢 Halo Orbits & Invariant Manifolds

Invariant manifolds about Halo orbits and their applications in the Circular Restricted Three-body Problem.

  • Joe Carpinelli

  • December, 2020

  • ENAE601 Final Project

Presentation Mode:
14.1 ms

📋 Project Overview

In the context of astrodynamics, manifolds are groups of trajectories that move toward or away from Lagrange points. To use invariant manifolds for interplanetary travel, several concepts must be developed and built on: lagrange points, periodic and quasi-periodic orbits within the Circular Restricted Three-body Problem, and manifolds about periodic orbits [1].

Outline

  • Brief review of Lagrange points

  • Periodic orbits (specifically, the subsection known as Halo orbits)

  • Finding Halo orbits (both analytically, and numerically)

  • Invariant manifolds about Halo orbits

Primary Reference

  • Megan Rund's Masters Thesis at California Polytechnic State University [1]

30.4 μs

Lagrange Points

  • Lagrange points are equilibrium points within the Circular Restricted Three-body Problem

21.0 μs

Stability at Lagrange Points

  • Like equilibrium points for all nonlinear systems, Lagrange points can be stable or unstable

  • We can find the stability of an equilibrium point by analysing eigenvalues of the Jacobian of the state vector

In the 2D case:

[ζ˙β˙ζ¨η¨]=[00100001UxxUxy02UyxUyy20][ζηζ˙η˙]

λ4+(4UxxUyy)λ2+UxxUyyUxy2=0

21.6 μs

Stability at Lagrange Points: Examples

  • Earth-Moon L4 is stable, and Earth-Moon L2 is unstable

)

593 ms

Periodic Orbits

Libration Orbit Families

  • Orbits about Lagrange points are known as Libration orbits

  • Lyapunov orbits are two-dimensional Libration orbits

  • Lissajous orbits are three-dimensional, semi-periodic Libration orbits

  • Halo orbits are three-dimensional, theoretically periodic Libration orbits

Finding Halo Orbits

  • The calculations required to solve for Halo orbits are numerically sensitive

  • We covered the procedure for iteratively solving for Halo orbits numerically in Lecture 16; this requires an initial guess

  • How do we find an initial guess for a Halo orbit?

38.2 μs

😇 Analytical Halo Solution

  • The Circular Restricted Three-body Dynamics can be derived using "Legendre poly- nomials", and there exists a third-order approximation (shown below) [5]

  • We can choose parameters to remove unstable (secular) terms from the expansion

Third Order CR3BP Expansion

Algorithm Inputs

  • Nondimensional mass parameter μ

  • Nondimensional Zaxis amplitude for the desired Halo orbit

  • Lagrange point to orbit (L1 or L2)

Algorithm Outputs

  • Position vector r0

  • Velocity vector v0

  • Estimated orbital period T

600 ms

Analytical Halo Examples: Earth-Moon L2

NOT numerically propagated!

798 μs

Analytical Halo Examples: Sun-Jupiter L1

NOT numerically propagated!

63.1 ms

🧮 Numerical Halo Solution

  • As discussed in Lecture 16, we can append the state transition matrix Φ(t0+ti,t0) to our state vector, and iteratively change initial conditions to numerically find a periodic orbit [1]

  1. Φ=I, x0=[x00z00y˙00Φ1Φ2Φ3Φ4Φ5Φ6]

  2. Propagate until y=0 again; Φ˙=FΦ where F=[0I3UXX2Ω], UXX is the matrix of second partial derivatives of potential U

  3. Calculate [δx0δy˙0]=([Φ41Φ45Φ61Φ65]1y˙[x¨z¨])1[x˙z˙]

  4. Set x0x0+δx0, y˙0y˙0+δy˙ and jump back to Step 1 (until x˙, z˙ are both within some tolerance of zero)

30.2 μs

🎢 Dynamics along Halo Orbits

Manifolds Exist

  • Each point along a Halo orbit is connected with an unstable manifold, and a stable manifold

  • The unstable manifold departs the Halo orbit, and the stable manifold arrives at the Halo orbit

  • How can we calculate the perturbation required to shift the spacecraft onto the manifold?

Finding Manifolds

  • We can use Eigenvectors of the Jacobian to calculate a state perturbation which will place the spacecraft onto a manifold at each point along the Halo orbit

  1. Propagate the Halo orbit for one period T, including the state transition matrix Φ(t0+ti,t0)

  2. Let the final state transition matrix be M; M=Φ(t0+T,t0)

  3. Calculate eigenvectors VSminreal(eig(M)), and VUmaxreal(eig(M))

  4. For each point i along the Halo orbit...

    • ViS=Φ(t0+ti,t0)VS, and ViU=Φ(t0+ti,t0)VU

    • XiS=Xi±ϵViS|ViS|, and XiU=Xi±ϵViU|ViU|

37.0 μs

Invariant Manifold Example

  • Pulled from Rund's Thesis [1]

  • Note that stable manifolds need to be propagated backward in time from the perturbation along the Halo orbit, becuase they converge on the Halo orbit [1]

143 ms

🪐 Manifold-based Transfer Designs

  • Design Overview summarized from [1]

Design Overview

  1. Find a desired Sun-Earth Halo orbit

  2. Place the spacecraft within the stable manifold of this Halo orbit

  3. Perturb the spacecraft from the Halo onto the unstable manifold to return towards Earth, and apply a maneuver to place the spacecraft on a Hyperbolic escape trajectory toward your destination planet

  4. At destination planet, apply a maneuver to place spacecraft onto stable manifold of destination Halo orbit

21.2 μs

🚀 Conclusions

Manifold Transfers

  • Lagrange points, and periodic orbits about Lagrange points are surrounded by collections of trajectories called manifolds

  • Manifolds can depart, or arrive at the Lagrange point / Libration orbit

  • Stable manifolds can bring us to Halo orbits for free (plus the cost to place the spacecraft onto the manifold)

Lessons Learned

  • Linear (Eigenvector) analysis works in Astrodynamics too

  • Calculations for iterating on Halo orbits are extremely numerically sensitive

  • Even with more complicated models (CR3BP), we need a patch-conic-like approach for interplanetary mission design

21.4 μs

📚 References

[1] Rund, M. S., “Interplanetary Transfer Trajectories Using the Invariant Manifolds of Halo Orbits,” , 2018.

[2] Howell, K. C., “Three-dimensional, periodic,‘halo’orbits,” Celestial mechanics, Vol. 32, No. 1, 1984, pp. 53–71.

[3] Vallado, D. A., Fundamentals of astrodynamics and applications, Vol. 12, Springer Science & Business Media, 2001.

[4] Richardson, D., “Analytical construction of periodic orbits about the collinear points of the Sun-Earth system.” asdy, 1980, p. 127.

[5] Lara, M., Russell, R., and Villac, B., “Classification of the distant stability regions at Europa,” Journal of Guidance, Control, and Dynamics, Vol. 30, No. 2, 2007, pp. 409–418.

[6] Koon, W. S., Lo, M. W., Marsden, J. E., and Ross, S. D., “Dynamical systems, the three-body problem and space mission design,” Free online Copy: Marsden Books, 2008.

[7] Williams, J., Lee, D. E., Whitley, R. J., Bokelmann, K. A., Davis, D. C., and Berry, C. F., “Targeting cislunar near rectilinear halo orbits for human space exploration,” 2017.

[8] Zimovan-Spreen, E. M., Howell, K. C., and Davis, D. C., “Near rectilinear halo orbits and nearby higher-period dynamical structures: orbital stability and resonance properties,” Celestial Mechanics and Dynamical Astronomy, Vol. 132, No. 5, 2020, pp. 1–25.

[9] NASA, NASA’s Lunar Exploration Program Overview, 2020.

[10] Carpinelli, J., “UnitfulAstrodynamics.jl,” https://juliahub.com/ui/Packages/UnitfulAstrodynamics/uJGLZ/, 2020.

13.2 μs

⌨️ Source Code

4.6 μs

Package Dependencies

  • The following packages were used – all are available in Julia's General package registry

11.4 μs
32.4 s

Finding Lagrange Points

6.0 μs
53.9 s

Unstable Lagrange Point (L2)

4.0 μs
17.2 s

Stable Lagrange Point (L4)

4.5 μs
2.0 s

Plot Analytical Halo Orbit

4.1 μs
386 μs

Northern Sun-Jupiter Halos

3.6 μs
7.3 s

Southern Sun-Jupiter Halos

3.8 μs
229 ms

Northern Earth-Moon Halos

3.9 μs
242 ms

Southern Earth-Moon Halos

3.7 μs
230 ms

Numerically Produced Halo

3.7 μs
61.7 s